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^ ' Abstract 

The ^^Be(p, d)^''Be reaction has been performed in inverse kinematics with a ra- 
dioactive ^^Be beam of E/A = 35.3 MeV. Angular distributions for the O"*" ground 
state, the 2+, 3.37 MeV state and the multiplet of states around 6 MeV in ^'^Be 
were measured at angles up to 16°^ by detecting the ^''Be in a dispersion-matched 
spectrometer and the coincident deuterons in a silicon array. Distorted wave and 
coupled-channels calculations have been performed to investigate the amount of 2+ 
core excitation in ^^Be^^. The use of "realistic" ^^Be wave functions is emphasised 
and bound state form factors have been obtained by solving the particle-vibration 
coupling equations. This calculation gives a dominant 2s component in the ^^Be^^ 
wave function with a 16% [2"'"0ld] core excitation admixture. Cross sections cal- 
culated with these form factors are in good agreement with the present data. The 
Separation Energy prescription for the bound state wave function also gives satis- 
factory fits to the data, but leads to a significantly larger [2"'"(8)ld] component in 
11 Be^,. 
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1 Introduction 



The nucleus ^^Be is of especial interest for several reasons. As is well-known, 
the ground state spin-parity is 1/2"*" in contradiction to the simple shell model 
and spherical Hartree-Fock prediction of l/2~. This "parity inversion" is cor- 
rectly predicted by, for example, recent psd-shell calculations of Brown [1]. 
The 2si/2 intruder orbital is lowered by the non-central part of the particle- 
hole interaction [2]. Moreover, ^^Be is often regarded as the classic one-neutron 
halo nucleus: the small single-neutron separation energy of 505 keV together 
with an assumed s-wave nature of the valence neutron leads to a very extended 
spatial distribution [3,4]. 

Several calculations of the ^^Be ground state structure have been performed. 
The theoretical approaches include: the shell model [5,6], the variational shell 
model [7], the Generator Coordinate model [8], and coupling of the neutron 
with a vibrational [9,10] or rotational core [11,12]. Most of these models cor- 
rectly reproduce the parity inversion and high-energy reaction data, but make 
very different predictions about the degree of coupling of an si/2 neutron to the 
^°Be 0+ ground-state core relative to a 6.^/2 neutron coupled to a 2+ excited 
core (the first excited state of ^°Be at 3.368 MeV). 

A direct test of the models for the structure of ^^Be^^ may be made by measur- 
ing the relative cross sections of one-neutron pick-up reactions feeding the O"*" 
and 2"*" states of ^°Be. Transfer cross sections depend on the overlap between 
the wave functions of the initial and final states through the radial neutron 
form factors uij{y). Standard distorted wave Born approximation (DWBA) 
analyses assume that these form factors are proportional to single particle 
wave functions U[J'{t), so that one may calculate cross sections independently 
of any prior assumption about the structure of initial and final states, apart 
from an overall normalisation factor. The latter is the spectroscopic factor, 
which is defined as the product of the overlap integral J ufj{r)r'^dr and a fac- 
tor (n + 1) [13], where n in the present case is the neutron occupation number 
of the 2sld shell in ^°Be. If one expresses the wave function of the 1/2"*" ^^Be 
ground state as the sum of the single particle and core excited components 

|"Beg, >= a|i°Be(0+) ® 2s > +/?|^°Be(2+) ® Id >, (1) 



the spectroscopic factors S(0"^) and S(2+) for transfer to the ground and first 
excited state of ^°Be should be directly related to a^ and /3^, respectively, 
assuming negligible population of the 2sld orbitals by ^°Be core neutrons.[^ 



^ Strictly speaking, a and (3 should be equal to the fractional parentage coefficients, 
the squares of which add up to unity. The relation between these and spectroscopic 
factors is given in the Appendix and Ref. [13]. 



0.55 


0.40 


0.42 


0.70 


- 


- 


0.78 


0.20 


0.20 


0.74 


0.19 


0.20 


0.80 


0.20 


0.20 


0.87 


0.10 


0.10 


0.92 


0.07 


0.07 



Table 1 

Coupling of the ^^Be ground state with the 0+ 0.0 MeV and 2+ 3.37 MeV states 
^''Be, as predicted by different structure calculations. The spectroscopic factors for 
the square of the wave function overlaps are given, together with the 2+ admixture 

S(0+® 2s) S(2+® Id) Rce 
Variational shell- model (Otsuka et al. [7]) 
PVM (Bhattacharya & Krishan [10]) 
CC rotational coupling (Nunes et al. [11]) 
Shell-model (Warburton & Brown [5,6]) 
Vibrational coupling (Vinh-Mau [9,14]) 
CC-DC (Esbensen et al. [12]) 
Generator-Coord. (Descouvemont [8]) 

Table 1 gives spectroscopic factors deduced from the various models cited 
above. These spectroscopic factors vary widely. For example, the standard 
Shell Model [5,6] predicts S(0+) = 0.74 and S(2+) = 0.19, while the Varia- 
tional Shell Model [7] gives S(0+) = 0.55 and S(2+) = 0.40. Also given in the 
table is the ratio Rce = S(2+)/[S(0+) + S(2+)], which provides a measure of 
the amount of the 2+ admixture in the wavefunction. This ratio shows large 
variations, although all models agree that the 2"*" admixture is smaller than the 
2s component. The dominance of the 2s component is also the conclusion of 
some recent experimental investigations [15,16]. However, the precise amount 
of core-excitation in the ^^Be ground state wavefunction remains uncertain. 

The present paper describes an approach to the ^^Be^s structure through the 
^H(^^Be, ^°Be)^H inverse-kinematics radioactive beam reaction. Data were also 
obtained for the ^H(^^N, ^^N)^H reaction at 584 MeV for comparison with the 
(p, d) reaction as measured in normal kinematics [17]. Our (^^N, ^^N) measure- 
ment provided an important cross-check of the analysis methods, particularly 
in regard to the calculation of the deuteron/heavy-ion coincidence efficiency 
and the ability to eliminate background from reactions on carbon. 

Brief accounts of the experiment have been presented in Refs. [18-20] . Prelimi- 
nary data were reported in Ref. [19] together with the ffist results of a standard 
DWBA analysis, which used single-particle form factors evaluated according 
to the usual Separation Energy (SE) prescription for the bound state wave 
functions. Spectroscopic factors deduced from that preliminary analysis [19] 
indicated a large [2+(8)ld] admixture in the ^^Be^^ wave function, exceeding 
most theoretical predictions. However, that interpretation was valid only to 
the extent that the radial wave function U£j{t) of the transferred neutron in 
^^Be, relative to a ^°Be core, may be approximated by the product of the SE 



single-particle form factor and a spectroscopic amplitude. This assumption is 
questionable in view of the large deformation parameter (/32=0.74 [21]) of the 
final nucleus ^"^Be, which could induce important coupling effects. These core- 
coupling effects may change the radial shape of the neutron wavefunction in a 
manner dependent on the angular momentum. This can directly affect the in- 
ferred spectroscopic factor deduced from transfer-reaction cross sections. The 
latter mainly depend on the squared amplitude in the surface region of the nu- 
cleus. In Ref. [20] we presented the results of DWBA calculations which used 
bound state form factors evaluated in the framework of the particle-vibration 
coupling model. These more sophisticated calculations indeed gave a signif- 
icantly smaller [2"'"(8)ld] admixture in the ^^Be^s wave function compared to 
the results with the SE wave functions. In the present paper, the vibrational 
coupling results are shown to be consistent with calculations which use radial 
form factors from the rotational excitation model of Nunes et al. [11]. In ad- 
dition, we present calculations which include coupling to inelastic channels, 
and discuss the most significant factors not included in the direct transfer 
calculations. 

The structure of the paper is as follows: After the experimental details (Sec- 
tion 2), the features of the spectra are discussed (Section 3). Next, we explain 
the extraction of the cross sections (Section 4), which is followed by a sec- 
tion on the analysis of the angular distributions with the distorted wave Born 
approximation (DWBA), including an estimation of the effect of coupling to 
inelastic channels (Section 5). The results are summarised and discussed in 
Section 6. 



2 Experimental Details 

2.1 Beams and targets 



The secondary ^^Be beam was produced by fragmentation of a 65 MeV/nucleon 
^^N beam from the GANIL cyclotrons, which bombarded a carbon target of 
thickness 1.03 g/cm^ located between the two superconducting solenoids of 
the SISSI device [22]. The ^^Be beam was analysed by the "alpha" spectrom- 
eter, which was operated as a fragment separator and had a 216 mg/cm^ 
aluminium achromatic degrader to reject unwanted ion species. After the de- 
grader, a small amount of He and Li isotopes remained in the beam which 
was 93% ^^Be. The magnetic rigidity of the second stage of the "alpha" spec- 
trometer and beam line after the degrader was 2.377 Tm, which corresponded 
to an average energy of 388.3 MeV {E/A = 35.3 MeV) for the ^^Be ions. 
With these settings, the intensity of the ^^Be particles was about 3 x 10"^ pps 
for 1.5/iA of ^^N, and the energy spread was 4.0 MeV. The full width at half- 



maximum (FWHM) spread of the incident beam angle at the reaction target 
was measured to be approximately 1.5° horizontally and 0.4° vertically 

The ^^Be beam was brought to a dispersive focus at the target, with a physical 
size of approximately 5 cm in width and 1.4 cm in height. By comparison, the 
aperture in the target frame was 10 cm by 3 cm. The principal target used 
in the experiment was a polypropylene (CH2CHCH3)„ foil of thickness 50 /im 
and density 0.896 g/cm^. The nominal target thicknesses were cross-checked by 
comparing the energy-loss of the ^^N direct beam (without degrading) with 
stopping-power and range tables [23]. Data were also collected with empty 
target frames during the course of the experiment, from which it is deduced 
that less than 2% of the counts in the ^°Be-gated focal plane spectra originated 
from the target frame or other sources, rather than from the polypropylene 
target. For the ^^Be beam, data on a 7.24 mg/cm^-thick carbon target were 
taken for background subtraction purposes in the singles spectra, as discussed 
in Section 3. 

For the ^H(^^N, ^^N)^H reaction, the ^^N primary beam was degraded in energy 
by a 745 mg/cm^ carbon foil in SISSI to 584 MeV (38.9 MeV/A), so as to 
more nearly simulate the properties of the ^^Be beam. The resulting spread 
in the incident beam angle at the reaction target was about 0.5° horizontally 
(FWHM) and 0.3° vertically (FWHM), which was somewhat smaller than that 
for the ^^Be beam. The size of the beam spot for the degraded-energy ^^N at 
the secondary reaction target was also smaller than that for the ^^Be, being 
approximately 4 mm wide and 4 mm tall at half maximum. 



2.2 Dispersion-matched ("energy-loss") magnetic spectrometer 



The energy spread of 4 MeV in the incident ^^Be beam, if not compensated 
in some way, would make it impossible to separate ^'^Be ejectiles in different 
excited states. The necessary compensation was achieved by the dispersion- 
matched magnetic spectrometer, SPEG [24], in which the initial 'analysis' 
beamline has a momentum dispersion = 9.86 cm/%. The spectrometer after 
the target is tuned to match the dispersed beam spot. In this way, all particles 
with the same reaction Q-value can be arranged to arrive at the same position 
along the focal plane. 

The SPEG spectrometer was placed at a central angle of 0°, with the ac- 
ceptance slits set to ±2.0° vertically and ±2.0° horizontally. The focal plane 
was instrumented to detect the ^°Be ejectile nuclei. Momentum and angle 
measurements were provided by two XY drift chambers [25] of the focal plane 
detector. This allowed the reconstruction of the position spectrum at the focal 
plane, which was at a slight angle to the normal to the mean particle trajec- 



tory. Unambiguous particle identification was achieved through a combination 
of the energy-loss signal from an ionisation chamber, the light output from a 
stopping plastic scintillator, and the time-of-flight measured by the scintillator 
relative to the cyclotron rf. 

The unreacted ^^Be beam was stopped in an active "finger" of plastic scin- 
tillator material placed at the high momentum side of the focal plane. The 
scintillator was coupled to a phototube and the light output and time signal 
were recorded in the acquisition system by a separate, downscaled trigger. 
This allowed the elimination of the He and Li contamination in the beam, as 
well as the electrons from the 13~ decay of the stopped ^^Be (half-life 13.8 s), 
and thus provided a true count of the number of incident ^^Be particles. 

The calibration of the focal plane in terms of the particle momentum (and 
hence excitation energy) was performed with a well-collimated ^^N beam inci- 
dent on a thin Au target with the spectrometer centered at 3°. The magnetic 
field of the spectrometer was changed in 1.5% steps, two above and two be- 
low the central field setting, and the centroid of the elastic scattering peak 
recorded. With the same beam and target, a mask with a pattern of holes was 
placed 65 cm after the target for the angle calibration, which was performed 
at each of the five field settings used for the momentum calibration. The hor- 
izontal and vertical angles at the focal plane, 6 foe and (pfoci respectively, were 
measured at each setting. The horizontal and vertical positions at the focal 
plane, Xfoc and Ufoo were also measured. With these calibration data, the hor- 
izontal angle at the target Otgt could be calculated directly from the measured 
9 foe- The optical properties of the spectrometer mean that 9 foe is independent 
of Xtgt and the other parameters. However, the vertical focussing is such that 
the vertical angle and the vertical position y are mutually dependent. This 
is an important effect for the present experiment, noting that (ptgt is equally 
as important as 9tgt in determining the reaction angle. Further, the vertical 
extent of the secondary beam spot at the target position, ytgt , is large for a 
secondary beam. To first order, the optical transport of SPEG in the vertical 
plane can be written as: 

Vfoc = (f) ytgt + (f ) (Ptgt 



[-) ytgt + [^) 4>tgt 



with the four coefficients depending on 9tgt and the momentum difference 5 = 
Ap/pQ measured relative to the central momentum, pq . These coefficients have 
been determined from the angle-mask calibration runs and the microchannel 
plate calibration described in Section 2.4. The measured yfoc and (pfoc values 
were then inverted to determine (ptgt and the less-important ytgt in the data 
analysis. The uncertainty in (ptgt was deduced by re-analysing the calibration 
runs with the final coefficients applied and was found to be a = 0.20° for both 



the ^^Be and ^^N runs, despite the smaller emittance of the latter beam. 
2.3 Beam tracking detectors 



The angles of the incident beam particles were measured event-by-event using 
two XY position sensitive drift chambers [26] located before the "analyser" 
magnet of the SPEG spectrometer, separated by 1 m from each other. This 
position was chosen in order to eliminate the products of any scattering or 
nuclear reactions inside the detector material. Each chamber consisted of four 
modules such that both X and Y were measured twice, but in opposite senses 
(to eliminate multiple coincident hits). The approximate overall dimension of 
each chamber was 7x7x7 cm^. They were filled with isobutane gas to a 
pressure of 20 mbar and each chamber had a total thickness of approximately 
0.65 mg/cm^ CH2 equivalent, including 4 /im of mylar in the gas- window and 
field-shaping foils. The timing signals from each module, relative to the focal 
plane plastic scintillator of the spectrometer, were fed into a multi-hit TDC. 
More than one X or Y signal from a given chamber could thus be recorded 
in a single event (triggered by the spectrometer focal plane) . Correlated pairs 
of X and Y times from each chamber were then selected by the data analysis 
software, by requiring that the sum of each pair should add up to a constant. 
If one measurement of an X-pair was missing in a given chamber, and if the 
time for the X-signal that was present was unique, the X-position could be 
recovered knowing that constant, although this was only required for a few 
percent of the total events. The combined efficiency of the two chambers for 
complete XY determination was approximately 88% for the ^^Be beam and 
96% for the ^^N beam. 

Since several magnetic elements in the beam line lay between the XY-tracking 
detectors and the target, the incident angles at the target needed to be cross- 
calibrated against the measurements at the tracking detectors. A tightly- 
collimated beam was sent through to the spectrometer focal plane with no 
reaction target in place. The currents in two horizontal and vertical steering 
magnets were then systematically changed, to steer the beam at different an- 
gles through the tracking detectors, and thence through the spectrometer. The 
difference signals {Xi — X2) and (Yi — Y2) from the tracking detectors were 
then related to the previously-calibrated Otgt and (ptgt from the focal plane 
detector. 



2.4 MicroChannel plate 



A position-sensitive microchannel plate detector [27] was placed approximately 
60 cm in front of the target and was used as a diagnostic device, as well as 



to calculate the aberration coefficients in Eqs. 2. It was not used for event- 
by-event correction of the deuteron angle because of poor detection efficiency. 
However, it was used to estimate the size of the beam spot, which in turn was 
used as input to the Monte-Carlo simulation program to calculate the coinci- 
dence efficiency as will be discussed in Section 4.1. The horizontal position was 
calibrated by stepping a tightly-collimated beam across the detector using the 
analysing magnet field. The vertical position was calibrated using a dispersed 
beam by illuminating a thick target which had a small hole in the centre, and 
stepping the target ladder up and down. A further calibration was performed 
with an a-source by placing a mask with a pattern of holes immediately in 
front of the detector target foil. 



2. 5 Detector array for coincident light particles 



The recoiling deuterons were detected in an array of ten position-sensitive 
sheet-resistive silicon detectors (CHARISSA) mounted in the target chamber. 
Each detector was 5x5 cm^ in area and 500 /im thick. These are the same 
detectors as used in the MEGHA array [28]. They were arranged around the 
space joining the extended beam spot and the spectrometer acceptance aper- 
ture, as shown in Fig. 1, and spanned angles between 5° and 35° relative to 
the central point on the target. The energy signals of the detectors were cali- 
brated with a 3-line a-source and the source peaks were monitored during the 
experiment. The position signals were calibrated with the a-source, by placing 
masks with a pattern of holes over the detector faces. The coincidence timing 
between the silicon detectors and the spectrometer focal-plane (taken from the 
plastic stopping detector) was set using the reaction -'^H(-'^^N, ^^N)^H ai E = 
584 MeV in kinematic coincidence, with the spectrometer moved temporarily 
to —2° and each CHARISSA detector moved in turn to +75°. 

The kinematics of the energy signal from the silicon detector as a function of 
the momentum of the ^°Be in the focal plane allows deuterons to be distin- 
guished from protons arising from the (p, d* ^ p + n) reaction. This is shown 
in Fig. 2. It was not possible to use a similar technique to separate the pro- 
tons from deuteron breakup in the (^^N, ^^N) data, because the higher beam 
energy in that case meant that the deuterons were not in general stopped in 
the silicon detectors. 

Being less sensitive to the incident beam angle [29], the deuteron angle proved 
useful for the beam angle calibration. Also, the vertical angle, after the ap- 
propriate kinematic transformation to the heavy-ion angle, could be used as 
an alternative to (ptgt discussed above. The further use of the deuteron angle 
was limited by the large size of the beam spot at the target. 
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Fig. 1. Schematic layout of the array of CHARISSA silicon detectors around the 
extended target. 

2.6 Electronics and data acquisition 



Standard NIM and CAMAC electronics were used to process the signals from 
the preamplifiers. The principal triggers for the standard GANIL VME-based 
acquisition system were: (i) a SPEG focal plane event, (ii) the downscaled 
beam in the plastic "finger", and (iii) a CHARISSA "singles" event. The latter 
was used only as a diagnostic. The data words for event-type (i) included a 
bit-pattern that indicated which, if any, CHARISSA detector had recorded a 
hit in coincidence. The time correlation was also recorded between the silicon 
signal and the heavy-ion in order to reject random coincidences. The data 
were recorded on magnetic tape and were independently analysed off-line at 
the IPN-Orsay and at the University of Surrey. 
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Fig. 2. Two dimensional plot of the energy deposited by particles in the CHARISSA 
array against ^''Be momentum in the focal plane. The upper groups correspond to 
deuterons from the (p, d) reaction. The more diffuse group below the deuterons 
corresponds to protons arising from the breakup of excited deuterons. The data 
include ^*^Be scattering angles from 0.4° to 1.2°. 

3 Discussion of Focal Plane Spectra 



3.1 ^'^N spectra 



Focal plane position spectra for the (-"^^N/^N) reaction at E/A = 38.9 MeV 
on the polypropylene target are shown in Fig. 3. The peak from the reaction 
on protons leading to the ^^N 1"*" g.s. dominates the spectrum. The experi- 
mental energy resolution is about 540 keV (FWHM). Calculations using the 
formulae given in Ref. [29] show that the combined uncertainty in the angle 
measurement and target effects contribute about 470 keV to this, leaving a 
contribution of about 270 keV from the spectrometer optics which corresponds 
to a momentum resolving power of approximately one part in 4300. 

The coincidence measurement (lower panel) removes the reactions arising from 
the carbon in the target, as well as possible background from ^^N — * ^^N + n 
breakup, seen in the singles spectrum (upper panel). Other states observed 
in i^N are the 0+ 2.313-MeV level and the second 1+ level at 3.948 MeV. All 
these states were seen in the "^^N(p,d) experiment of Ref. [17], and with a 
similar ratio of strength to the present data. In the singles spectrum (upper 
panel), weak yield to low-lying states in ^^C from the ^^C(^^N,^''N)^^C reaction 
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Fig. 3. ^^N focal plane spectra: Upper panel in singles, and lower panel in coincidence 
with deuterons in the CHARISSA array. The spectra are for ^^N laboratory angles 
from 0.3° to 0.6°. 

is observed. The peak immediately to the left of ^^C ground state at about 
3.9 MeV excitation, and of similar width to the ground state peak, is assigned 
to ^^C excited to 3.854 MeV rather than ^^N excited to 3.948 MeV. The 
latter assignment is unlikely, since the peak would be broadened by inflight 
7-emission to an estimated width of about 1.36 MeV. The excitation energy 
of ^'^C (possibly mutual with ^^N) for the broad peak near channel 225 is 
~ 9 MeV. 



3.2 ^^Be spectra 



The ^°Be focal plane position spectra for reactions with the ^^Be beam were 
accumulated at two different magnetic field settings separated by 0.9%, for a 
total number of 3.7 x 10^ incident ^^Be particles. Only the data for the first 
field setting, amounting to roughly 6/7*^ of the total, were used in the final 
analysis. 
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Fig. 4. ^''Be focal plane spectra. Upper panel: in singles; Lower panel: in coincidence 
with deuterons in the CHARISSA array. The spectra are for ^'^Be laboratory angles 
from 0.4° to 1.2°. Superimposed on the upper panel is a spectrum taken on a carbon 
target, normalised to the same number of beam particles and equivalent carbon 
thickness as that for the polypropylene target. 

The (^^Be, ^°Be) singles spectrum is shown in the upper panel of Fig. 4. The 
low-lying states of ^^C from ""^^C in the target are at most weakly observed; 
as a reference point, the ground state for ^^C is expected at channel 222. 
In contrast, there is a strong yield to states at high excitation, such as the 
peak at channel 165, corresponding to an excitation in ^^C of ~ 10 MeV, or, 
more likely, to mutual excitation of ^^C to ~ 3.5 MeV and of ^"^Be to 6 MeV. 
This strong yield contaminates the peaks from the reaction on hydrogen. The 
carbon-scattering origin of these counts is proved both by their absence in 
the spectrum taken in coincidence with deuterons in the CHARISSA array 
(Fig. 4, lower panel) and their presence in the singles data taken on a pure 
carbon target (Fig. 4, overlay on upper panel). The preferential population of 
high excitation states in -"^^C (with an assumed 5/2^ structure) agrees with Q- 
value and angular- momentum matching considerations [30]. The coincidence 
spectrum (Figure 4, lower panel) shows clean separation from other counts for 
the ground state and the first excited state (3.368 MeV, J'^ = 2+) of ^"^Be. 
As observed in the case of the ^^N beam spectra, the coincidence measure- 
ment very effectively removes the background which principally arises from 
the stripping reaction on the carbon in the target. 

The energy resolution for the ground state peak is about 700 keV. The angle 
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uncertainty and target straggling accounts for only about 250 keV of this, 
from which we deduce that the contribution from the spectrometer optics was 
some 650 keV. This corresponds to a momentum resolving power of approxi- 
mately one part in 1200, which is considerably worse than that observed for 
the (^^N, ^^N) reaction and is presumably a consequence of the larger beam 
spot size in the (^^Be, ^'^Be) case. 

All the peaks corresponding to ^°Be excited states in the spectra are broadened 
by gamma decay of the nuclei in flight. The measured fwhm of the peak for the 
3.37 MeV state is 1.14 ± 0.29 MeV. The maximum spread (i.e. base width) 
from the recoil-broadening is 2f3E^ = 1.89 MeV, where jS = v/c = 0.28 is 
the average ^'^Be velocity. If one assumes an isotropic angular distribution 
of 7-rays, giving a sin(^) dependence to the spread, the estimated fwhm is 
1.26 MeV, in agreement with observation. Full treatments of the effect on the 
broadening when the angular distribution of the 7-emission is known are given 
by Beene and Devries [31] and Pelte and Schwalm [32]. 

The dominant peak in the ^°Be spectra arises from transfer to a group of four 
closely-spaced levels near 6 MeV excitation. Kinematic matching [30] does not 
especially favour the 6 MeV region in ^'^Be, thus the strong yield to this region 
reflects a strong overlap of initial and final states. One explanation of the 
intense "6- MeV peak" is the pickup of one P3/2 neutron in ^°Be core, feeding 
the (si/2,pjw2) 1~ and 2~ states at 5.960 and 6.263 MeV. There might also 
be smaller contributions of neutron pickup to the second 2"*" and O"*" states 
at 5.958 and 6.179 MeV. However, an alternative, although not necessarily 
incompatible explanation invokes the two-centre shell-model [33,34]. In the 
two-centre shell-model, the ^^Be nucleus is described as two alpha particles 
and three orbiting neutrons [35]. The removal of a neutron should leave ^"^Be in 
a similar two-centre orbit configuration. In the case of ^°Be the most strongly- 
coupled configurations are not associated with the ground state but with the 
states around 6 MeV [36,37]. We note also that the two-centre shell model 
can also account for the inversion of the s- and p- states resulting in the 
1/2+ ground state in ^^Be. This is because of the lowering of the Q'^ = 1/2+ 
component of the deformed d5/2 orbitals. 

In the singles spectrum, the peak at 6 MeV is superimposed on a background 
corresponding to the high energy tails of the (^^Be^^^Be + n) and (d— ^p + n) 
break-up reactions near thresholds. The residual nucleus ^°Be is unbound to 
particle emission above 6.812 MeV, so no further peaks are expected above 
this multiplet. As previously discussed in Section 2.5, the contribution of the 
(p, pn) reaction to the ^°Be spectra has been removed by energy conditions in 
CHARISSA which select coincident recoil deuterons. 
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4 Extraction of Cross Sections 



In this section, we present a description of the simulation program used to 
estimate the particle detection efficiencies, which is needed for the extraction 
of the cross sections from the deuteron coincidence data. All the cross sections 
presented in this paper are those from the analysis of the coincidence spectra. 
Despite the larger number of counts in the peaks, the cross section extraction 
from the singles spectra has no better statistical accuracy than that from the 
coincidence spectra, because of the need to subtract a background which is not 
accurately determined. Analyses of some singles spectra have been performed 
as a cross-check of the estimation of the coincidence efficiency. 



4-1 Detection efficiency 



Data in coincidence with the CHARISSA detectors were corrected for the vari- 
ation of geometrical deuteron detection efficiency as a function of the heavy-ion 
(^°Be or ^^N) laboratory angle. This efficiency was calculated by the use of an 
extended version of a Monte-Carlo simulation program [38], which accounts 
for the emittance and beam spot size of the secondary beam, and for the un- 
certainty of the scattering angle determination. Two efficiencies need to be 
taken into account in the calculation of the coincidence differential cross sec- 
tion. Firstly, the efficiency of detecting deuterons in the CHARISSA array for 
a given heavy-ion detection angle. Secondly, since the beam spot is large, (par- 
ticularly for the ^^Be beam), not all the heavy ions at an angle less than the 
nominal spectrometer slits setting (which assumes a point beam spot) will be 
transmitted. There is a further complication arising from the large variation in 
incident beam angle. Thus, a second efficiency factor needs to be calculated, 
which is the ratio of the number of heavy-ion events generated at a given angle 
divided by the number of counts accumulated in the corresponding "detected" 
angle bin. The angular resolution of the detection system is also taken into 
account in the program. This latter is an important effect at small angles, and 
partly explains the apparently high deuteron detection efficiency near 0°. 

Figure 5 shows the calculated efficiencies for the ^^Be beam. The beam spot 
size and envelope were set to simulate the experimentally-measured values 
given in Section 2.1. 

The procedure for the extraction of the coincidence cross sections is to in- 
tegrate the cleanly-separated peaks and divide by the total efficiency calcu- 
lated over the appropriate angular range. For several forward angles in the 
(^^Be,^°Be) data, where the background from the carbon target is relatively 
small, it was possible to extract the ground state yield in the singles spectra. 
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Fig. 5. Detection efficiency calculated by the Monte-Carlo program of deuterons in 
the CHARISSA array (dotted line) and of ^^^Be^s through the spectrometer slits 
(dashed line) as a function of ^'^Be laboratory angle. The solid line is the combined 
overall coincidence efficiency. Calculations at discrete angles are plotted, connected 
by straight lines. 

The ratio of the coincident to singles yield for the ^°Be ions measures the 
deuteron deuteron efficiency and can be compared to the one calculated by 
the simulation program. For an angular range between 2.5 to 8.5°„, where the 
calculated efficiency curve is roughly constant with angle (Fig. 5) the aver- 
age coincident to singles yield ratio was found to be 0.28 ± 0.04. This agrees 
within errors with the average efficiency of 31.4% calculated by the simulation 
program. 



4-2 N Angular Distributions 



Cross sectional angular distributions for the 1^ ground state of ^^N are shown 
in Fig. 6. These are from the analysis of the deuteron-coincident spectra. Other 
states in ^''N were not analysed because of the low yield and the uncertainty 
in subtracting the underlying background. 

The last angle bin extracted is for 6iab = 2.1 — 2.3°, corresponding to a mean 
angle of 6cm ~ 29°. In principle, one could go beyond the nominal angular 
acceptance limit of the spectrometer (±2° square for a point beam spot, giv- 
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Fig. 6. Experimental angular distributions for the present inverse kinematics re- 
action to the ground state in ^^N (open circles) compared to the (p, d) data of 
Snelgrove and Kashy [17] (triangles). The curve is to guide the eye only. 

ing a diagonal limit of 2.83°) because of the spread in incident beam angle. 
However, the fall-off in reaction cross section and lack of coincidence efficiency 
at large angles prevents this. 

In Fig. 6 we also show the (p,d) data of Snelgrove and Kashy [17], which 
were measured at Ep = 39.8 MeV. Our inverse-kinematic data at E/A = 
38.9 MeV are in reasonable agreement with the (p, d) data. We conclude that 
our analysis technique in general, and our simulations of the deuteron efficiency 
in particular, are reliable at least up to a centre-of-mass angle of 24° for 
p^i5jyj^ ^^N)d. For that reaction, 24°^ corresponds to a deuteron laboratory 
angle of approximately 28°. If we take the latter as a limit of reliability on 
the deuteron scattering angle, and apply it to the p(^^Be, ^^Begs)d reaction, 
we obtain a rough limit of 15°„ up to which we are confident in the efficiency 
simulation. The corresponding limits for the reaction leading to the 3.368 MeV 
and 6.1 MeV excited states of ^"^Be are 18°^ and 21°^, respectively. 



4-3 Be Angular Distributions 



The final experimental angular distributions for the 0^ and 2+ states and the 
6 MeV multiplet are shown in Fig. 7 and tabulated in Table 2. The error 
bars are statistical. For the four most-forward angle points, the systematic 
errors arising from uncertainties on the angle determination in SPEC and 
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Table 2 

Differential angular distributions for ^^Be(p, d) leading to states in ^''Be. 



0+ 


0.00 MeV 


2+ 3.34 MeV 


6 MeV states 


"cm 


a{e) 


"cm 


a{e) 


"cm 


a{e) 


(deg) 


(mb/sr) 


(deg) 


(mb/sr) 


(deg) 


(mb/sr) 


2.50 


5.21 ±1.39 


2.62 


4.85 ± 1.30 


2.73 


19.54 ±2.54 


5.08 


5.16 ±0.73 


5.39 


6.03 ± 0.74 


5.62 


20.85 ±1.32 


8.27 


4.05 ±0.53 


8.66 


4.87 ±0.55 


9.00 


17.80 ± 0.98 


11.60 


1.30 ±0.27 


12.15 


4.74 ± 0.49 


12.63 


22.08 ±1.04 


15.02 


0.55 ±0.25 


15.74 


4.57 ±0.79 


16.39 


19.47 ±1.42 



in the CHARISSA detectors, and on the variation efficiency of the detection 
system as function of the scattering angle, are estimated to be of the order 
of 15%. The largest angle points have been revised from the data shown in 
ref. [20] as a result of a new estimation of the coincidence efficiency [39]. These 
largest-angle points are close to the limits of the reliability criterion deduced 
from the (""^^N, ^^N) data above, and are more sensitive to the approximations 
made for the beam profile, giving an estimated 25 — 30% additional systematic 
uncertainty from the efficiency calculation. The largest-angle points were not 
used in the extraction of spectroscopic factors in ref. [20] , neither are they so 
used in the present paper. 



5 Analysis of Angular Distributions 

5. 1 Optical model potentials 



Different combinations of optical potentials for the entrance and exit channels 
have been tried in the calculations presented below, in order to test the sensi- 
tivity of the extracted spectroscopic factors to the input parameters. All the 
optical potentials used in the present analysis have the standard Woods-Saxon 
or Woods-Saxon derivative form. 

For the entrance channel, three principal optical potentials have been used. 
The most-recent global nucleon-nucleus optical parameterisation is the "CH89" 
one of Varner et al. [40]. This has dependences on energy, mass and isospin, 
adjusted for a range of stable nuclei from masses A = 40 to 209. However, data 
from recent proton elastic scattering experiments involving ^°Be and ^^Be ra- 
dioactive beams (at similar incident energies to the present work) were only 
reproduced with the CH89 parameterisation if the depth of the real Saxon- 



Table 3 

Optical model parameters for the p + ^^Be entrance (P„) and d + ^'^Be exit (Dn) 
channels. Potentials D2 and -D3 are adiabatic potentials taken from the nucleon 
potentials of Ref. [40] and [43] , respectively, for use in the Johnson-Soper approxi- 
mation. Potential depths are in MeV while diffuseness and reduced radii are in fm. 
Radii are defined in terms of the reduced radii as Rx = r-x x A^'"^. 



Label Ref V r„ a„ W„ W^ r,,, a„, V, 



so ' so 



Pi 


[40] 


41.0 


P2 


[42] 


56.0 


P3 


[43] 


58.4 


Di 


[44] 


76.8 


D2 


a) 


73.8 


D2 


b) 


74.8 


D3 


a) 


106.: 


D3 


b) 


107.: 



1.15 0.69 3.64 7.77 1.14 0.69 5.9 0.80 0.63 

1.08 0.64 3.70 3.91 1.25 0.68 5.6 1.01 0.60 

1.12 0.57 - 10.39 1.12 0.50 5.5 1.12 0.57 

1.15 0.81 - 23.04 1.34 0.68 - 

1.15 0.72 4.08 15.83 1.14 0.72 5.9 0.78 0.63 

1.15 0.72 3.71 16.17 1.14 0.72 5.8 0.78 0.63 

106.3 1.13 0.60 - 15.9 1.13 0.53 5.5 1.13 0.57 

107.3 1.14 0.60 - 16.1 1.14 0.53 5.5 1.13 0.57 

a) adiabatic potential for d + ^'^Be^^. 

b) adiabatic potential for d + ^°Be2+. 

Woods well were reduced by a factor of 0.88 [41]. The CH89-derived potential 
Pi in Table 3 corresponds to p + ^^Be at 35.3 MeV with the real well depth ad- 
justed by this factor. Secondly, the 1980 parameterisation of Fabric! et al. [42] 
was used to generate the proton potential P2, assuming a deformation param- 
eter for ^^Be of zero. Thirdly, the proton-nucleus potential P3 was obtained 
from the parameterisation of Watson, Singh and Segel [43] , which was derived 
from the analysis of elastic scattering data in the Ip shell.0 

The situation for suitable exit channel scattering potentials is not as satis- 
factory as that for the entrance channel. A "global" parameterisation may 
be found in Ref. [44], but it has been derived from the analysis of deuteron 
scattering on a limited number of A > 24 targets. No polarised beam data 
were available at the time that the analysis was made, and so a spin-orbit 
potential is lacking. This potential, adapted for 36 MeV deuterons, is listed 
as Di in Table 3. Zwieglinski et al. [45] also used the Perey and Perey [44] 
parameterisation in the analysis of their ^''Be(d, p)^^Be data but added their 
own spin-orbit potential (depth 7 MeV). Test calculations that we have per- 
formed with such an added ad-hoc spin-orbit potential to Di gave only small 
differences to the cross sections in the present, rather restricted, angular range 
of interest. 



'^ We note that the spin-orbit parameterisation in Ref. [43] is non-standard, having 
l/Rso df/dr instead of the usual 1/r d//dr. Since the spin-orbit potential has little 
effect on the angular distributions, we have nevertheless kept to the standard form. 
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5.2 Single-particle form factors for a Woods- Saxon well 



Theoretical differential cross sections <Jdw have been calculated using the zero- 
range DWBA code DWUCK4 [46]. Test calculations were made with correc- 
tions for the effects of finite range and the non-locality of optical potentials; 
in these tests, a finite range parameter of 0.621 and non-locality ranges of 
0.85 and 0.54, for nucleons and deuterons, respectively, were used. Differences 
in the cross sections with and without the finite range correction were found 
to be less than 10%. The effect of the non-locality correction on the spectro- 
scopic factors was found to be generally less than 20%. For some particular 
optical potentials, the cross sections changed by about 30%, although in these 
cases the effect on the O"*" and 2"*" cross sections is correlated to some extent. 
The DWUCK4 calculations presented here have no non-locality corrections. 
Transitions to the 0+ ground state and 2+ first excited state were assumed to 
proceed through neutron pick-up in the 2si/2 and ld5/2 orbitals, respectively. 
The doublet of 1~ and 2~ states around 6 MeV was assumed to be excited 
through a P3/2 neutron pick-up from the ^'^Be core of ^^Be. 

For the calculations described in this section, bound state neutron form factors 
were calculated in a Woods-Saxon well according to the usual Separation 
Energy (SE) prescription [13]. This consists of adjusting the well depth so 
that the eigenvalue of the Schrodinger equation is equal to the experimental 
separation energy, thus ensuring correct asymptotic behaviour of the wave 
function. In the calculations a spin-orbit Thomas term was used with A=25 and 
two different geometries for the Woods-Saxon-well: (1) standard geometrical 
parameters ro=1.25 fm, a=0.65 fm, and (2) ro=1.15 fm, a=0.57 fm, from the 
parameterisation of Ref. [43] for p-shell nuclei. 

Different sets of optical potentials were used in the generation of the distorted 
waves in the entrance and exit channels, as detailed in Section 5.1. Attempts 
to increase the radius parameter of P2 and P3 by 25% to account for the 
particularly large matter radius [47] of ^^Be were not successful in reproducing 
the present angular distributions, and the geometry given in Refs. [42,43] was 
adopted for the production run calculations. 

Two fundamentally-different types of deuteron potentials have been used. Po- 
tentials deduced from elastic scattering analyses, such as the Perey potential 
Di [44], are suitable for standard DWBA calculations. On the other hand, 
the adiabatic deuteron breakup approximation (ADBA) proposed by Johnson 
and Soper [48] is generally known to improve the description of (p,d) and 
(d,p) reactions, by accounting for the effects arising from the break-up of the 
deuteron in the nuclear field. Such adiabatic deuteron potentials were obtained 
by folding the neutron and proton global potentials at half the deuteron en- 
ergy, according to the prescription of Ref. [49] . Potential D2 was derived from 
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the CH89 [40] parameterisation and potential D-^ from the Watson, Singh and 
Segel [43] parameterisation. 

These different combinations of optical parameter sets reproduce reasonably 
well the rapid decrease of the ground state /=0 cross section with angle and 
the rather flat angular distribution observed for the 2+ state and the 6 MeV 
peak (see examples in Fig. 7). The data in general were poorly reproduced 
with the global potential [50] extrapolated from (d, d) scattering on A > 27 
nuclei, or the potentials used for analysing ^'^Be(d, d) data at 12 and 15 MeV 
[21]. (The results from these "failed" calculations are neither shown, nor are 
the potentials considered further.) The calculated angular distributions for a 
given Ij transfer were normalised to the four most forward angle data points 
by a least square fit procedure, in order to determine the spectroscopic factors 
S. These factors were deduced from the relation a^xp = NSa^w/i'^J + l); 
where j is the total spin transfer, with the (p, d) zero range normalisation 
factor A^ of 2.29 [46]. The square of the isospin Clebsch-Gordan, C^, which is 
sometimes included in the definition of S, is unity for ^^Be^s ^^°Be + n. 

The spectroscopic factors 8(0'*') and S(2+) are shown in Fig. 8 for various sets 
of optical potential parameters. The experimental S(2+) values obtained by 
assuming a pure 6.^/2 neutron pick-up change by less than a few per cent if 
instead an / = 2 neutron pick-up in the higher-lying 6.3/2 subshell is considered. 
The absolute spectroscopic factors strongly depend on optical parameters, 
with for example a factor of about 1.8 between the S(2+) values obtained 
with the Percy and Percy deuteron potential Di and the adiabatic potential 
D3. The sum of spectroscopic factors for 0^ and 2'^ states, related to the 
occupation number of the s-d shell in ^^Be^s, was found to vary between the 
extreme values of 0.74 and 1.80, with an average value of 1.09. On the other 
hand, the ratio Rce = 8(2"'")/ [8(0"'") +8(2"'")] related to the admixture of core 
excited components in the wave function, is less dependent on the calculation 
parameters, as shown in Fig. 8. 

The average value oi Rce deduced from this 8E analysis is 0.51, with a standard 
deviation of 0.09. Note that the lowest values of Rce are obtained from the 
ADBA calculations with the optical parameter set P3-D3. These parameters 
also provide the best fit to the data for angles below 14°^, and when 8E form 
factors are used with either geometry (1) or (2) for the binding potential, 
one obtains Rce values of 0.30 and 0.32, respectively. We therefore deduce a 
conservative lower limit of 30% core excitation admixture in the ^^Be^^ wave 
function from the analysis which used single-particle form factors calculated 
in a Woods-8axon well with the 8eparation Energy procedure. 

The angular distribution for the 6 MeV peak is well reproduced by calculations 
for / = 1 transfer to the 1~ and 2~ states at 5.96 and 6.26 MeV (see Fig. 7), 
without consideration of a possible couphng of ^^Bcgs to the 23" and O2" states 
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Fig. 7. Theoretical angular distributions calculated under the DWBA obtained with 
single-particle SE form factors for states in ^'^Be. The points are the experimental 
angular distributions. 

unresolved in the multiplet. With the optical model parameter set P3-D3 and 
form factor geometry (1), the summed spectroscopic factor is found to be 1.40. 
This summed spectroscopic factor is in good agreement with the results of the 
shell model calculations of Warburton and Brown [5,6] in an extended lp-2s- 
Id basis, which predict spectroscopic factors of 0.69 and 0.58 for calculated 1~ 
and 2~ states at 5.96 and 6.23 MeV. However, we recall that this parameter set. 
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Fig. 8. Spectroscopic factors and the ratio i?ce = S(2^)/[S(0^)+S(2^)], extracted 
from a standard DWBA analysis, using different optical parameter sets and sin- 
gle-particle form factors SE-1 and SE-2, calculated with geometrical parameters (1) 
and (2), respectively. 



P3D3, gave the lowest spectroscopic factors for the transition to the ^'^Be 3.37- 
MeV state (see fig. 8). Calculations for the 1~ and 2~ states with parameter 
set P3-D1 gave a summed spectroscopic factor of 2.69, which is twice the shell 
model predictions. On the other hand, the shape of the experimental angular 
distribution suggests that the amount of any / = or / = 2 contribution 
is quite small. The agreement of the P3-D3 calculations with the shell model 
predictions is thus circumstantial evidence supporting that choice of optical 
model potentials. 
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5.3 Effects of core coupling on transfer cross sections 



The comparison in the previous section between experimental and theoretical 
spectroscopic factors is based on two assumptions: (i) a single-particle bound 
state form factor, and (ii) negligible interference from two-step processes. We 
have noted in the Introduction the particularly large experimental deforma- 
tion parameter /52=0.74 [21], of the 2"*" state in ^°Be as deduced from (p,p') 
data. This might induce significant coupled channel effects on the (p,d) cross 
sections. In addition, the SE assumption of the proportionality of the bound 
state form factors to the single-particle wave functions may be a poor approx- 
imation if certain basic conditions are not fulfilled. These conditions include 
having an experimental binding energy close to that expected from the sim- 
ple shell model and negligible residual interaction between the core and the 
transferred nucleon [13]. In the case of the pick-up of the halo neutron orbiting 
around a ^°Be core, coupling to the excited core may also modify the radial 
shape of the neutron wave function in a significant way. 

A description of the theoretical formalism and of the calculation of the vi- 
brational form factors is given in the Appendix (further details will be given 
in Ref. [39]). Briefly, the ^^Be^^ wave function was written in the form given 
in eq. (1), and a set of coupled equations were solved for the nucleon wave 
functions. Calculations were performed for both sets of Woods-Saxon well ge- 
ometries (1) and (2) previously adopted for the SE calculations, and for two 
different conditions on the deformation parameter value j32- By adjusting the 
strength of the real central and spin-orbit potential, the eigenvalues of the 
coupled equations were made equal to the experimental separation energies. 

For the VIB-1 and VIB-2 cases listed in Table 4, the experimental f32R value 
of 1.84 fm from the analysis of ^'^Be(p,p') data [21] was used. The values 
of the potential depths are given in Table 4. Note the unusually large Ko 
(about twice the standard spin-orbit value) needed to reproduce the ^^Bej,^ 
parity inversion. A similar result was found by Nunes et al. [11]. In spite of 
the different radial shapes of uij{r), similar theoretical spectroscopic factors 
gt/i ^ g g^ ^^ gt/i ^ Q j^g deduced from both calculations. 

Sl/2 «5/2 

A smaller core-excitation (S^'^ ~ 0.10), is predicted by the calculations VIB- 
Ib and VIB-2b, which have a more standard, smaller spin-orbit strength, cor- 
responding to a Thomas-Fermi A factor of 25. This requires an increased de- 
formation parameter (32 of approximately 0.9 to reproduce the experimental 
binding energies. 
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0.84 


0.67 


0.16 


0.17 


0.89 


0.65 


0.09 


0.13 




(0.66) 




(0.28) 


0.83 


0.80 


0.16 


0.26 


0.87 


0.79 


0.11 


0.21 




(0.79) 




(0.38) 



Table 4 

Parameters and resulting spectroscopic factors of vibrational coupling calculations 
for ^^Be. Spj ,2 measures the overlap between ^^^Be^s and l/2~ state at 0.32 MeV. 
Ss-^,,2 and Sdg.j are for transitions to the 0^ ground state and first 2^ state of 
^'^Be, respectively. S^^^ values were obtained from the present DWBA analysis with 
optical parameter set PsD^, using the normalised Uij wave functions as form factors. 
Numbers given in parentheses correspond to the results obtained using standard SE 
form factors in the same well geometry. Potential depths are in MeV and distances 

in fm. 

^2 yp Vso rp a S|,^^2 S*4^2 Sf^P S*^^^^ S"^^^^ 

VIB-1 0.68 51.5 15.4 1.25 0.65 0.84 0.84 

VIB-lb 0.92 51.0 7.0 1.25 0.65 0.61 

SE-1 1.25 0.65 

VIB-2 0.74 60.3 13.3 1.15 0.57 0.81 

VIB-2b 0.90 60.0 8.3 1.15 0.57 0.67 

SE-2 1.15 0.57 

5.3.1 Direct one-step transfer 

Theoretical cross sections for direct transfer using external form factors were 
calculated with the code DWUCK4 [46]. The use of vibrational form factors 
(normalised to unity) in place of the SE ones in DWBA or ADBA calcula- 
tions neither appreciably modifies the shape of angular distributions nor the 
magnitude of the 0^ cross section, but enhances the theoretical cross section 
for the 2"*" state by typically a factor of two. As shown in Ref. [20], the effect 
of the coupling on the radial wave functions is only minor for Si/2, but signifi- 
cantly shifts the 6.5/2 form factor outwards, increasing its magnitude by about 
40% in the important asymptotic region outside the nucleus. This is respon- 
sible for the strong difference between the cross sections calculated with the 
vibrational coupling model and the SE method. The omission of these cou- 
pling effects must cast some doubt on the physical meaning of the S(2"'") and 
Rce values deduced from the SE method, as plotted in Fig. 8. Experimental 
spectroscopic factors Sfj^^ determined as in Section 5.2 but using vibrational 
formfactors are given in Table 4. The corresponding R^e ratios, giving the 
experimental amount of core excitation, range between 0.17 and 0.24, to be 
compared with the minimum value of 0.30 from the standard SE analysis. The 
S(0'*") and S(2"'") values obtained using SE form factors (Section 5.2) are also 
given (within parentheses) in Table 4 for comparison. 

Alternatively, one can use the unnormalized form factors as ingredients in 
the DWBA calculations and directly compare the predicted cross sections to 
experimental data. This is done in Fig. 9, which shows the results obtained 
with optical parameter set P3-D3, and vibrational form factors corresponding to 
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different well parameters, without any further renormalisation to the data (i.e. 
the Sjj factors are implicitly used in these calculations). As was the case for 
the SE method in Section 5.2, the calculated cross sections agree well with the 
revised last angle data point in the 2"*" distribution. Note that the calculated 
2+ cross sections are not proportional to the theoretical spectroscopic factors, 
since the VIB-1 and VIB-2 cross sections differ by a factor 1.5 for the same 
S^Jl value (0.16). The best agreement between experimental and theoretical 
cross sections for the 2"*" state is obtained using form factors VIB-1, calculated 
with the large spin-orbit parameter (Vso = 15.4 MeV) needed to reproduce 
the energy of the excited l/2~ state with the same Woods-Saxon well. By 
contrast, the VIB-lb and VIB-2b results (with standard spin-orbit strengths 
and [32 values of ~ 0.9) underestimate the 2+ cross section by a factors of 1.4 
and 1.9, respectively. 

Direct transfer calculations were also performed with radial form factors from 
the core excitation model (CEM) of Nunes et al. [11]. This model assumes 
a rotational structure for the ^'^Be core and possible population of the 2si/2, 
ld5/2 and ld3/2 by the valence neutron. We have taken spectroscopic factors 
of 0.85, 0.13 and 0.02 for these levels, respectively (parameter set Bel2-b from 
Ref. [52]). The predicted cross sections, again with the optical parameter set 
P3D3, are shown as dashed lines in Fig. 9. They are rather similar to the 
results of the vibrational form factors (VIB-1) and reproduce the data fairly 
well. 

The Rce value deduced from the present analysis of the p(^^Be, ^'^Be)d cross 
sections therefore provides evidence for a dominant (~ 80%) coupling of the 
halo 1/2+ ground state of ^^Be to the 0+ ground state of ^"^Be. The absolute 
gs -^ gs spectroscopic factor depends on the geometry adopted for the neutron 
well, as might be expected, but it was found not to be substantially modified 
by the channel coupling treatment of form factors (unlike the case of the 2"*" 
transition). Considering the uncertainties from optical potentials alone, the 
present (p, d) experimental value (0.67 or 0.80, for geometrical parameters 1 
or 2, respectively), may be considered in good agreement with the theoretical 
predictions S*'' of the vibrational coupling model in Table 4. 

Both of our experimental and theoretical results for the ground state transition 
are also in good agreement with the results of two previous (d, p) experiments 
which used radioactive ^°Be targets. These experiments gave spectroscopic 
factors of 0.73 (Ref. [21]) and 0.77 (Ref. [45]) for incident deuteron energies 
of 12 MeV and 25 MeV, respectively. One can also compare the (d, p) values 
for the transition between the ^^Be^s and the first l/2~ state at 0.32 MeV 
in ^^Be, 0.63 (Ref. [21]) and 0.96 (Ref. [45]), with our calculated values S^';^^ 
which range between 0.61 and 0.84 (Table 4). 
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Fig. 9. Theoretical cross sections for single step transfer with vibrational coupling 
form factors calculated with different geometrical and spin-orbit parameters. Results 
obtained with form factors from the Core Excitation Model (CEM) of Ref. [11] are 
also shown. The calculations are not normalised to the data. 

5.3.2 Contribution of two- step processes to the 2^ cross section 

An important question is whether processes involving inelastic excitation in 
^°Be and ^^Be could contribute to modify the 2"*" cross section in a significant 
way, as did the coupling in the form factor. This has been checked by cal- 
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Fig. 10. Coupling schemes for calculations involving inelastic excitation. Top panel: 
scheme for the ZR one-way coupling calculations leading to the 2+ state in ^^Be. 
Bottom panel: scheme for the full FR-CRC calculations discussed in Section 5.3.3 

culating two-step amplitudes for inelastic excitation followed by transfer (or 
vice- versa) with the zero range code CczR [53], which are then added to the 
direct transfer amplitude (with the appropriate phases). The coupling scheme 
considered in these calculations is shown at the top of in Fig. 10. 

Only one-way coupling was considered for inelastic excitation, as the inelastic 
coupling strength /32R was taken equal to the experimental value in ^*^Be [21] 
extracted with this assumption. Vibrational coupling form factors VIB-1 have 
been used for the Si/2 and d5/2 neutron transfers. In addition to the two-step 
(p, dd') excitation of the 2"*" state via the 0"*" ground state of ^°Be, possible 
two-step (pp',d) paths via 3/2"'" and 5/2"*" states in -^-^Be with [2"'"(8)Si/2] con- 
figuration, were also considered. These 3/2"*" and 5/2"*" states are predicted to 
have an excitation energy similar to that of the ^°Be 2"*" state (3.37 MeV) by 
the weak coupling limit of the vibrational model. Experimentally, good candi- 
dates for these weak coupling partners of the 2"*" state are the unbound levels 
at 2.69 and 3.41 MeV, with possible J'^ values 3/2+ or 5/2"^ [54]. Inelastic 
coupling strengths and Si/2 neutron form factors of the (p, dd') path were also 
adopted for (pp',d), with the appropriate phases and amplitudes. Results are 
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Fig. 11. Cross sections for the 2^ state in ^''Be calculated with CczR. "direct" is for 
the transfer route alone (no coupling), "(pp',d)" and "(p, dd')" are for the individual 
two-step routes (including inelastic one-way coupling), "ZR-DW2" is for both the 
two-step routes together with the direct. Results for vibrational form factor VIB-1 
and two different optical parameter sets are shown. 

shown in Fig. 11 for optical parameter sets P3-D3 and P3-D1, which gave a good 
description of the O"*" differential cross sections calculated with direct transfer. 

For calculations involving inelastic scattering explicitly, where the macroscopic 
model of a deformed scattering potential is used for the inelastic excitation 
form factor, the standard d + ^°Be optical potential Di is suitable. On the 
other hand, adiabatic deuteron potentials are not intended for this application, 
as these potentials a priori reproduce neither the relevant elastic nor inelastic 
scattering channels. The results labelled P3-D3 in Fig. 11 have therefore been 
obtained with the adiabatic deuteron optical potential D^ as the generator of 
the distorted waves in the (d+^'^Be) channel, whilst the inelastic form factors 
were proportional to dV /dr, the derivative of the standard deuteron optical 
potential Di [44]. 
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The enhancement of the 2"*" cross section at forward angles from the two-step 
processes considered here, strongly depends on the optical parameter sets used 
in the calculations, as they induce large differences in both the direct transfer 
and the inelastic scattering amplitudes. For P3-D3, interferences with two-step 
processes modify the slope of the angular distribution, but do not change 
the cross section at forward angle by more than 15% relative to single step 
transfer, so that it is still in agreement with the present data. For P3-D1, it 
is seen in Fig. 11 (top panel) that interferences with the (pp', d) and (p, dd') 
paths induce a large enhancement of the 2"*" cross section (by a factor of 2.4 
at 0°) relative to direct transfer. 



5.3.3 Comparison of full with one-way coupling calculations and effect on the 
0^ cross section 

Further coupled-reaction-channels (CRC) calculations have been performed 
with the finite-range program FRESCO [55]. These were intended partly as a 
test of the approximation of restricting the coupling to one-way, as done in 
the calculations of the previous section. In addition, since the de-excitation 
of the 2"*" level is now considered, the effect of dynamical coupling on the 
0^ cross section can be examined. Finite-range calculations, however, are not 
compatible with the Johnson-Soper ADBA approach which takes into account 
deuteron breakup effects, since the adiabatic potential used is obtained from 
the three-body (core-proton-neutron) Schrodinger equation assuming zero sep- 
aration between the proton and neutron. Hence the calculations presented in 
this section use the standard scattering optical potentials P3-D1. For simplicity, 
the (pp'd) path through the hypothetical 5/2+ state of ^^Be was not included 
in these Fresco calculations (that through the 3/2'^ state was considered). 
Lastly, all the calculations in this section included the "remnant term" in the 
interaction potential, normally assumed to cancel in zero-range distorted-wave 
programs (the effect of the inclusion of the remnant term is discussed below). 
Thus, for several reasons, the calculated two-step angular distributions are 
not expected to be similar to those presented in the previous section that 
used the same optical model potential combination. The couplings considered 
are shown schematically at the bottom of Fig. 10 

In Fig. 12 Fresco calculations are presented for full CRC, for one-way cou- 
pling and for the direct transfer alone. The bottom panel of the figure shows 
that for the transfer leading to the 2"*" state in ^°Be, the difference between the 
one-way coupling (FR-DW2) and the full coupling (FR-CRC) cross sections is 
no more than 10%. As previously seen in the ZR calculations of Section 5.3.2, 
there is a large difference between the one-way coupling and the direct transfer 
cross sections for this optical model choice. 

As regards the transition to ^°Be O'^ ground state, the FR direct, one-way 
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Fig. 12. Cross sections for the 0"'' ground and 2"*" 3.37 MeV states in "'^''Be calculated 
with Fresco. "FR-CRC" is for the fuh finite-range CRC, "FR-DW2" is for one-way 
inelastic coupling together with direct transfer, "direct" is for the transfer alone (no 
coupling). 



coupling (now including the 2^ to O"*" de-excitation) and full coupled channels 
angular distributions all show obvious differences in magnitudes, with the one- 
way coupling distribution having in addition a small angular displacement in 
the minimum compared to the others (Fig. 12, top panel). At small angles, 
the one-way coupling cross section is about 60% of the direct cross section. 
The FR-CRC cross section magnitude is approximately midway between the 
two others. At large angles (> 20°), the CRC cross section becomes the lowest 
of the three. 

One should be aware of a possible problem with these comparisons, in that the 
same optical potentials have been used for the direct and coupled-channels cal- 
culations rather than refitting the parameters (and the /?2 coupling strength) 
to elastic and inelastic scattering data with the same model, had such data 
existed. Nevertheless, it seems clear that coupling to inelastic channels could 
in principle affect not only the ^°Be 2"*" and 0^ cross sections, but, more im- 
portantly, the ratio Rce- 
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5.3.4 Effect of the remnant terms in the interaction potential 

As mentioned in Section 5.3.3 above, one difference between the finite-range 
Fresco and zero-range Dwuck programs, is that the former allows the in- 
clusion of the full complex remnant terms (also known as the "indirect" part) 
in the interaction potential. To be explicit, the prior representation of the 
interaction for a pickup reaction 5(6, a) A with a = h + x is [56] 

W^ = VaA -Ua = V,A + {VbA - t/a), (3) 



where V^n is an effective nucleon-nucleus interaction and Ua is the distorting 
potential used for the entrance channel wave Xa- The binding interaction V^a 
is the only part from Eq. 3 retained in DwuCK, it being argued that there is 
considerable cancellation between the "remnant" terms (Ha --f^Q)-[3 However, 
the cancellation can never be complete, and the omission of these terms negates 
the post-prior equivalence of the interaction. 

From a comparison of the results of one-step Fresco calculations, the cross 
section for the O"*" ground state increases by about 20% when the remnant 
terms are included whereas that of the 2"*" excited state is not significantly 
changed. Thus the ratio of the 2+ and 0+ cross sections is affected, without 
the radial wave functions. The remnant terms appear to take into account the 
radial extent of the 0+ wave function being greater than that of the 2"*". 



5.4 Effect of ^^Be recoil excitation and breakup 



A new method for analyzing (p,d) reactions involving a weakly bound halo 
nucleus (C-|-n) has been developed recently by Timofeyuk and Johnson [57] 
as an extension to the ADBA. Their adiabatic approach treats the effect of 
the recoil excitation and break-up (REB) of the halo nucleus on the transfer 
cross section, by replacing the distorted wave in the p+(C+n) channel by 
an effective distorted wave, derived from the potential (p-l-C) according to 
a definite prescription. The distorted wave function in the (d+C) channel is 
obtained as in the ADBA, to account for the effect of deuteron break-up on 
neutron transfer. 

The application of this method to the ^^Be(p,d)^°Be reaction at 35 MeV shows 
that REB enhances forward angle cross sections of the 0+ and 2+ states, lead- 
ing to an appreciable reduction of corresponding spectroscopic factors (mul- 
tiplied by about 0.54 and 0.71, respectively) [57,58]. Such reductions applied 



^ In fact the zero-range approximation would make the calculation of these terms 
difficult since each is a function of different coordinates. 
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to our results would give a slight increase of Rce and would not change the 
lower limit determined for the amount of core excitation in ^^Be^s- The anal- 
ysis method used in Refs. [57,58] leads to an R^e value close to 0.5, which falls 
within the range of values shown in fig. 8. However, it has to be observed that 
the absolute value of the spectroscopic factor for the 0"*" state obtained by Tim- 
ofeyuk and Johnson even without the REB effects, S(0+) = 0.35 ± 0.04 [58], 
is only in agreement with the lowest of our SE analyses, namely, that with 
the optical potential P1-D3 (note that refs. [57,58] used the SE-1 geometry 
for the radial form factors). Their absolute value for 8(0"'") of 0.19 ± 0.02 with 
REB is difficult to reconcile with, e.g., the measurement of total reaction cross 
section [59] and that for high-energy one-neutron removal [16]. 



6 Discussion and Conclusions 



The ^H(^^Be,^°Be)^H reaction has been investigated in order to provide insight 
on the structure of the ground state of the halo nucleus ^^Be. The analysis 
has attempted to relate the large cross section observed experimentally for the 
excitation of the 2"*" state at 3.37 MeV with the amount of [2"'"(g>ld] core-excited 
component in the ground state wave function. 

An important ingredient in the calculations is the choice of method for calcu- 
lating the radial wave function of the transferred nucleon in the ^^Be nucleus. 
A common choice, which we have labelled SE, is to calculate the wave func- 
tion in a Woods-Saxon well, with the depth adjusted to reproduce the known 
separation energy. Another important ingredient is the use of adiabatic exit 
channel optical potentials to approximate the effect of the deuteron breakup. 
We find that any (p,d) reaction analysis using single particle form factors in 
the SE method will result in a ^'^Be 2"^ core excitation admixture >30%. This 
result would be in agreement with the predictions of variational shell model 
(42%) [7] although not with most other theoretical models. 



However, the validity of the SE single-particle form factors has been brought 
into question by coupled-channel calculations in the framework of the particle- 
vibration coupling model. The radial wave function of the (\^/2 transferred 
neutron is found to be strongly modified by the interaction with the ^'^Be 
deformed core, enhancing the cross sections relative to the classical SE predic- 
tions. One-step transfer cross sections calculated with form factors from the 
present vibrational coupling approach, or from a core excitation model based 
on the assumption of a rotational ^"^Be core [11], are in agreement with the 
data, and predict only ~10 - 20% admixture of the [2'^(8>ld] configuration in 
^^Be^s. Such a dominance of the s-wave component is in agreement with many 
theoretical calculations (e.g., Refs. [5,11,9,12]). 
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Because of the large deformation of ^''Be, a valid concern is that coupling 
channel effects may also contribute to the transfer cross sections in a signifi- 
cant way, via two-step (inelastic + transfer) processes. This has been tested 
by calculations in Sections 5.3.2 and 5.3.3 using rather simple assumptions. 
The effect of the inelastic coupling could indeed be significant, but it is depen- 
dent on the optical model potentials. The breakup of the deuteron seems to 
have a much more dramatic effect, even to the extent of changing the conclu- 
sions that may be drawn, regarding the importance of two-step processes. This 
strongly indicates that more complete calculations, including a more elaborate 
treatment of breakup, are required. These are beyond the scope of this work. 
Effects of ^^Be recoil and break-up should ideally also be included, as recent 
calculations [57] have shown that they could affect transfer cross sections sig- 
nificantly. We have also shown a sensitivity of the calculations to the inclusion 
of the remnant term in the interaction potential. This appears to affect only 
the 0+ cross section (at the 20% level), and consequently the extracted ratio 
of spectroscopic factors, Rce is also modified. The necessity of performing both 
elastic and inelastic scattering experiments in the (p -|- ^^Be) and (d -|- ^''Be) 
channels at the appropriate incident energy is evident. These are needed to fix 
the parameters external to the transfer process. It should also be noted that 
it is the nucleon-core potential, p -|- ^"^Be, which is required as input to the 
ADBA and REB adiabatic calculations. 

Three other key experimental results concerning the structure of ^^Be^s should 
be considered. Following a demonstration by Suzuki et al. [60] of the sensitiv- 
ity of the ground state magnetic moment of ^^Be to the [2+® ld5/2] component 
of the wave function, a measurement of the magnetic moment has been made 
at ISOLDE and has given a value of — 1.6816(8)/iN [15]. This result implies 
a relatively pure 2si/2 halo state with essentially no ld5/2 admixture, but a 
precise interpretation is dependent on the assumed quenching of the single par- 
ticle magnetic moment. Another recent experiment was a high-energy single- 
neutron removal ("knockout") reaction [16] which used a ^^Be beam. The 
authors' measured cross sections to the particle-bound levels in ^°Be agreed 
well with single-particle cross sections calculated in an eikonal model [61] mul- 
tiplied by spectroscopic factors from the shell model [5] . It was concluded that 
^^Be ground state has a dominant 2s single particle character with a small 
Id component. The third experiment is the measurement of the total reaction 
cross section of ^^Be [62], which was one of the earliest indicators of the ex- 
tended halo structure of this nucleus. Reaction calculations [63] can interpret 
this in terms of the average radius of the halo, which in turn is related to the 
percentage of s-wave component in the ^^Be wave function. The result [59] is 
that the s-wave fraction is of order 0.8 to 0.9, and larger d-wave admixtures 
could not be reconciled with the measured cross section. 

In conclusion, we have explored the structure of ^^Be ground state via a single- 
neutron transfer reaction with a ^^Be radioactive beam. From the experimental 
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point of view, the data presented here constitute a successful implementation 
of the technique of transfer reactions in inverse kinematics to study the struc- 
ture of light exotic nuclei, as suggested in, e.g., Refs. [29,64-66]. The present 
"best estimate" of the ^^Be^g wave function is a dominant 2s component with a 
0.16 [2^® Id] core-excitation admixture. This value may be better defined by 
future calculations that principally incorporate a model of deuteron breakup 
(as well as that of ^^Be) within a CRC framework. 
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Appendix A: Transfer form factors in a vibrational coupling ap- 
proach 



It has long been established that in the presence of core-polarisation admix- 
tures created by an interaction Hamiltonian Hint-, radial form factors of one- 
nucleon transfer reactions should no more be approximated by the product of 
a spectroscopic amplitude and a single-particle wave function U^^ deduced 
from the standard Separation Energy procedure, but have to be determined 
by solving the Pinkston-Satchler coupled equations [67]. In the present anal- 
ysis of the ^^Be(p, d)^°Be reaction, these coupled equations were solved in the 
framework of the particle- vibration coupling model [68,69], with the program 
CCVIB [51]. The interaction Hamiltonian was chosen as: 

H^nt = -R^ E ^=T (4m + (-1)"^^C,_m) r,V(^, 0) (.1) 



where c^^ and Cl-m are operators for phonon creation and annihilation and 
Pl is the deformation parameter. Here the term dV/dr corresponds to the 
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deformation of the central part V{r) of the potential, with no action on the 
spin-orbit part. Coupled-channel calculations were performed for both 1/2"'" 
ground state and 1/2"" first excited state at 0.32 MeV, with the additional con- 
straint that the same well depth Vq reproduced the experimental separation 
energies relative to the ground state of ^°Be, thus ensuring the right asymp- 
totic behaviour of the wave functions. This reproduction of the separation 
energies could be obtained by adjusting either the value of the (32 parameter 
or the strength of the spin-orbit part of the potential. The configuration space 
was truncated to the Ol, and 2f core states in ^'^Be (first order vibrational 
couphng), coupled with one neutron in the lp3/2 or lpi/2, and 2si/2 or ld5/2 
orbitals, for the 1/2" and 1/2"*" states, respectively. 

The resulting wave functions were normalised for each final state according 
to the prescription J{J2ijufAr))r'^dr = 1. The relative weights of the differ- 
ent configurations are given by the corresponding squared amplitudes afj = 
{J ufj{r)r'^drY. Strictly speaking, spectroscopic factors Sij are defined [13] as 
the product of the squared overlap integral {J Uij'uij{r)r^dry and {n + 1), 
where n is the occupancy of the orbital Ij in the core nucleus and t/f^ is the 
corresponding single-particle wave function calculated in the same Woods- 
Saxon potential well. The factors Sf'^ given in Table 4, which are in fact the 
squares of parentage amplitudes although we call them theoretical spectro- 
scopic factors for convenience, were calculated by assuming n = and, for the 
particular cases of Si/2 and d5/2, by approximating the unbound single-particle 
wave functions [//?', by the bound f/^'^^ radial form factors determined by the 
SE method. They were found to differ from the relative weights af, by less 
than 1%. 

The radial wave functions Uij{t) obtained by solving coupled equations can be 
directly used in the calculation of transfer cross sections by DWBA and/or 
CRC methods, so that the comparison of theoretical and experimental re- 
sults provides a test of the validity of the present vibrational coupling model 
and other ingredients such as optical potential parameters. Another possi- 
ble procedure is to normalise the calculated wave functions by dividing them 
by /■u?(r)r^(ir, and use the renormalised Uij{t) as "realistic" form factors in 
DWBA analysis, thus giving experimental spectroscopic factors expected to 
be more reliable than those extracted by the conventional SE method. Both 
methods were used in the present analysis of the ^^Be(p,d)^''Be reaction. 
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